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1. INTRODUCTION 

We shall demonstrate that the Fermi surface of dense neutron matter may expe- 
rience a rearrangement near the onset of pion condensation, due to strong momentum 
dependence of the effective interaction induced by spin-isospin fluctuations. In par- 
ticular, the Fermi surface may take the form of a partially hollow sphere having a 
spherical hole in its center. Thus, a second (inner) Fermi surface may form as high- 
momentum single-particle states are filled and low-momentum states are vacated. 
The influence of this phenomenon on the superfluid transition temperature of the 
Fermi system is characterized with the help of a separation transformation [1,2] of 
the BCS gap equation. 

The present effort can be viewed as a revival of the search for physical realizations 
of abnormal occupation in infinite, homogeneous Fermi systems - plausible instances 
in which the quasiparticle distribution differs from that of an ideal Fermi gas and 
Fermi liquid theory breaks down. This search (broadened to Bose systems) was in 
full swing at the time of the Third Workshop on Condensed Matter Theories, held in 
Buenos Aires in 1979, and it was led by Workshop founders Valdir Aguilera-Navarro, 
Ruben Barrera, Manuel de Llano, and Angel Plastino, with important contributions 
from James Vary, John Zabolitzky, and others [3-5]. 



2. THE BUBBLE REARRANGEMENT 

Rearrangement of the characteristic Landau quasiparticle distribution signals a 
breakdown of Fermi-liquid theory and can lead to profound alteration of the ortho- 
dox behaviors we have come to expect based on this pervasive physical picture. For 
the sake of simplicity, our considerations will be limited to homogeneous systems, for 
which the Fermi-liquid distribution uf(p) coincides with the momentum distribution 
of an ideal Fermi gas. The actual quasiparticle distribution n(p) must deviate from 
n F(p) = 9(pf ~ p) if the necessary condition for its stability is violated. At zero 
temperature, this condition states that the change of the ground-state energy Eq re- 



mains positive for any admissible variation 5n(p) away from n F (p). More specifically, 
stability of a given quasiparticle distribution implies 

SE = Jt(p,n(p))5n(p)-^>0, (1) 

where £(p,n(p)) = e (p,n(p)) — p is the quasiparticle energy measured with respect 
to the chemical potential p. For n(p) = n F (p), the condition (1) is violated if e(p) 
rises above p, at p < p F , or if e(p) falls below at p > pp. A rearrangement of 
quasiparticle occupancies is precipitated when the density p attains a critical value 
p c F at which the relation 

Z(p,n F (p);p cF ) = (2) 

exhibits a bifurcation leading to a new root p = po- This relation usually serves only 
to determine the Fermi momentum pp. 

In homogeneous systems, the simplest type of rearrangement of the momentum 
distribution n(p) of quasiparticles of given spin and isospin maintains the property 
that its values are restricted to and 1, but the Fermi sea becomes multiply con- 
nected (cf. Refs. [4-7]). In particular, we may suppose that at densities exceeding 
the critical value p cF , the normal-state distribution 9(p F — p) is altered by the for- 
mation of a "bubble," or particle void, over a range pi < p < pi < pp, with the 
Fermi momentum pp readjusted to maintain the prescribed density. This distri- 
bution is formally represented by n(p) = 9(pi — p) + 9(p F — p)8(p — pi), where as 
usual 9(x) = 1 if x > and vanishes otherwise. One then has three Fermi surfaces: 
two inner surfaces located at pi and pi, along with the usual outer surface at pp. 
However, a more dramatic rearrangement can also occur, resulting in a distribution 
with partial occupation of quasiparticle states that lacks the distinctive trademark 
of Fermi-liquid theory, namely a discontinuity of n(p) at the Fermi surface. In this 
scenario, called fermion condensation, there exists a finite range of momenta over 
which the quasiparticle energy coincides with the chemical potential, corresponding 
to the creation of a "fermion condensate" [8-12]. 

Any change of n(p) from the normal-state distribution n F (p) must entail an 
increase of the kinetic energy of the quasiparticle system. Accordingly, the anticipated 
rearrangement only becomes possible if it is accompanied by a counterbalancing 
reduction of potential energy, which implies that the effective interaction between 
quasiparticles has acquired a substantial momentum dependence. The emergence of 
such a strong momentum dependence is exactly what one expects to occur as the 
density p is increased toward the critical value p c for a second-order phase transition 
in which a branch of the spectrum u s (k) of collective excitations of the Fermi system 
collapses at a nonzero value k c of the wave vector k. 

To justify this expectation, we follow Dyugaev [13] and consider the behavior 
of the quasiparticle scattering amplitude -F(pi, P2, k) = z 2 T(p ll p 2 ; k, oj = 0)M*/M 
near of the phase-transition point. Here r(pi, p 2 ; k, uj) is the ordinary (in-medium) 
scattering amplitude, M* is the effective mass, and z is the renormalization factor 
specifying the weight of the quasiparticle pole. The amplitude F can be written as 
the sum F r + F s of a regular part F r and a singular part F s , with the latter taking 
the universal form 

Ks;p y (Pi, P2, k; p -> p c ) = -O aS O^D(k) + O a yOf, S D (|pi - p 2 + k|) (3) 



in terms of the propagator D(k) of the collective excitation. This form has been 
derived with due attention to the antisymmetry of the two-particle wave function 
under exchange of the particle coordinates (spatial, spin, isospin). The collective 
propagator is conveniently parametrized according [13] 



where the parameter (3(p), with (3{p c ) = 0, measures the proximity to the phase- 
transition point. The vertex O appearing in Eq. (3) determines the structure of the 
collective-mode operator and is normalized by Tr(OO^) = 1. Specifically, the choice 
O = 1 is made in treating the rearrangement of the quasiparticle distribution due to 
collapse of density oscillations [14], while O = a is appropriate when studying the 
rearrangement of npip) triggered by the softening of the spin collective mode [15]. 
In the present investigation we will be concerned with dense, homogeneous neutron 
matter in which abnormal occupation is induced by spin-isospin fluctuations; thus 
the pertinent operator is O = (a ■ k)r. 

Cutting through the details, the most important features of the model defined 
by Eqs. (3) and (4) are that the function F s (pi, p 2 , k = 0) ~ -D(pi — P2) depends on 
the difference pi — p 2 and that in the neighborhood of the soft-mode phase-transition 
point this dependence becomes very strong. 

Eqs. (3) and (4) furnish a suitable basis for efficient evaluation of the single- 
particle spectrum £(p) near the second-order phase transition. We exploit a straight- 
forward connection between and the scattering amplitude F(pi,p 2 ,k = 0), 
thereby circumventing the awkward frequency integration that would be encountered 
in an RPA approach. This connection is made through the relation 



which may be derived by means of the Landau-Pitaevskii identities [16-18]. The 
contribution to Eq. (5) from the singular part (3) of F can be easily integrated over 
the momentum p to obtain 



In stating this result, we assume that the contribution to £(p) from the regular, 
nonsingular part of F can be simulated by replacing of the bare mass M appearing 
in Eq. (5) by an effective mass M*. The generally accepted values for this effective 
mass are in the range 0.7-0.8 for the pertinent densities in the neutron-star interior. 

According to Migdal and collaborators [19,20] (see also Ref. [21]), a dramatic 
phase transition can occur when the density p of neutron matter in the liquid core of a 
neutron star reaches a critical density p C7r of some 2-3 times the equilibrium density 
Po of symmetrical nuclear matter. The spin-isospin collective mode collapses at a 
finite wave vector k = k c ~ pp and a phase transition identified as pion condensation 
sets in. A prominent feature of the ground state of the system beyond the phase- 
transition point is the presence of a condensate of spin-isospin density waves. It 



D-\k)=p 2 +1 \k 2 /kl-l)\ 



(4) 







should be clear from Eqs. (3) and (4) that spin-isospin fluctuations with k ~ k c 
are strongly amplified in the neighborhood of the transition as a consequence of the 
divergence of the propagator D(k — > k c , p c ). 

We are thus led to apply Eq. (6) to dense neutron matter close to the onset of 
neutral pion condensation. Employing the parametrization (4) in Eq. (6), we obtain 
the working formula 

p 2 1 f 1 d^pi 

^ = 2M* + 2~J f3 2 + 7 2 ((p - Pl ) 2 - k^/ki n{Pl) (2^)3 ' (7) 

Ideally, one would like to extract quantitatively reliable values for the input param- 
eters (3, 7, and k c from a convincing ab initio treatment of neutron-star matter. 
Unfortunately, no such treatment is yet available. Moreover, current predictions of 
the critical density p c * range from 0.2 to 0.5 fm -3 (i.e., some 1-3 times po), depending 
on what theoretical assumptions are implemented [19-21]. 

This situation leaves us with no alternative but to carry out calculations for 
several representative choices the parameters of the microscopic model. Inserting the 
formula (7) into Eq. (2), we determine the critical density p cF at which the solution 
of Eq. (2) bifurcates. For p > p cF this equation then yields two new momenta pi and 
pi at which vanishes, and between which £(p) is positive. The bubble region lies 
between these two momenta. 

The necessity for brevity in this presentation precludes the explicit presentation 
of the numerical results that have been derived for the spectrum and for the 
phase diagram of dense neutron matter. The neutron spectrum has been calculated 
at the critical densities p c F corresponding to three different sets of model parameters: 
(a) 7 = 1.25m*, k c = 0.9 p F , (3 2 = 0.22m 2 (p cF ~ 1.19 po), (b) 7 = 1.25m*, k c = 
0.9 p F , [3 2 = 0.25 m 2 (p cF ~ 1.76 /9 ), and (c) 7 = 1.25m*, k c = p F , f3 2 = 0.13m 2 
(pcF — 1-88 po), where m* is the pion mass. The choice k c = 0.9p F for the critical 
wave number is suggested by earlier numerical investigations [20]. Two different 
positions were found for the bifurcation point, namely po = (for parameter sets (a) 
and (b)) and p ~ 0.12 p F (for set (c)). The phase diagram of neutron matter in 
the variables p (measured in po) and [3 2 (measured in m 2 ) has been constructed at 
k c = 0.9p F for four different values of 7 (1.0, 1.2, 1.4, and 1.6, in units of the pion 
mass). Plots of the results of these calculations may be found in Refs. [18,22]. 

Variation of the parameters /3, 7, and k c within sensible bounds can significantly 
affect the phase diagram and hence the extent, in density, of the phase with rear- 
ranged quasiparticle occupation. Even so, our numerical study has revealed four 
characteristic and generic features of the bubble rearrangement. 

(i) The critical density p cF for the rearrangement is less than the critical density 
p c * for pion condensation. Since both phenomena arise from the strong mo- 
mentum dependence of the amplitude F(p 1? p 2 ;k — > 0), rearrangement of the 
quasiparticle distribution may be regarded as a precursor of pion condensation. 

(ii) The bifurcation point associated with formation of a bubble in the neutron mo- 
mentum distribution is located at small momenta, po < 0.2p F , regardless of the 
applicable value of p c *. 

(iii) The spectrum £(p) exhibits a deep depression for p ~ (0.5 — 0.6)p F . 



(iv) The ratios Pcf/Pc-k and Pq/pf are insensitive to the actual value taken by p C7T 
within the range of plausible theoretical predictions. 

The emergence in neutron matter of one or more new Fermi surfaces positioned 
at low momentum values would provide a new avenue for rapid direct-Urea neutrino 
cooling of neutron stars [22] . More broadly, the creation of new Fermi surfaces by the 
mechanism we have described - as well as the more profound rearrangement involved 
in fermion condensation - would call for revision of many of the conclusions that have 
been developed within Fermi liquid theory. Here we shall focus on some elementary 
properties of pairing in the reconfigured system. 



3. PAIRING IN THE PRESENCE OF ABNORMAL OCCUPATION 

We assume, for the sake of simplicity, that beyond the instability point there exist 
only two Fermi surfaces, an outer one corresponding to the usual Fermi momentum 
Pf and an additional inner one at pi, lying close to the origin in momentum space. 
Hence we consider the limiting case pi = in our original specification of the "bubble" 
rearrangement. Also for the sake of simplicity, we restrict the analysis to 1 So pairing, 
for which the BCS gap equation has the familiar form 

A(p) = - J V(p,p 1 )£- 1 (p;T)A(p 1 )dT ll (8) 

where V is the effective particle-particle interaction and we employ notations dr = 
p 2 dp/2w 2 for the volume element and S~ 1 (p;T) = [2E{p)]~ l tanh [E(p)/2T] for the 
usual combination of tanh temperature factor and energy denominator 2E(p). The 
appearance of the gap function A(p) in the superfluid quasiparticle energy E(p) — 

1/2 

[£ 2 (p) + A 2 (p)] renders the gap problem nonlinear. The quantity £(p) is to be 
interpreted as the single-particle spectrum in the system with pairing turned off. 

Adopting the strategy for solving gap equations that was introduced in Ref. [1] 
and elaborated in Refs. [2,3], we write the block V, identically, as a separable part 
plus a remainder that automatically vanishes on the outer Fermi surface. Hence we 
write 

V(pi,P2) = V f <Pf(pi)(Pf(P2) + W(p 1 ,p 2 ) (9) 

and take 4>f(p) = V(PiPf)/Vf, where Vf = V(j>f,Pf)- It follows directly that 
W(p,pf) = W(pf,p) = 0, as required. In the ordinary case where there is only one 
Fermi surface, this decomposition allows us to replace the singular nonlinear inte- 
gral equation (8) by two equivalent equations: (i) a nonsingular quasilinear integral 
equation for a T-independent shape factor x(p) = A(p)/Ap- and (ii) a nonlinear 'al- 
gebraic' equation for the T-dependent gap value Ap(T) = A(pp,T). In the present 
case where there are two Fermi surfaces, we must extend the procedure of Ref. [1] 
to deal consistently with the inner Fermi surface as well as the outer one. This is 
accomplished by decomposing the interaction term W appearing in Eq. (9) in the 
same manner as before, setting 



W{p u p 2 ) = W/0/(pi)0/(p 2 ) + Y(pi,p2) 



(10) 



with 0j(p) = W(p,pj)/W/ and W> = W(p/,pj) = V(pj,pj) - V 2 (>f,P/)/V(>f,Pf), 
so that = Y(pi,p) = = pi?) = 0. The above relations entail 

the boundary values 

0f(PfO = 1, 0/(pj) = 1, 0j(pf) = O, (11) 

while the key quantity 4>f(pi) oc V(pj,p_f) describes the connection between the 
quasiparticles of the two Fermi surfaces in the particle-particle channel. If 4>f(pi) 
vanishes, these surfaces are disconnected and the problem is trivialized. 

In the general case where V(pi,pf) 7^ 0, substitution of Eqs. (9) and (10) into 
the BCS gap equation (8) gives 

A(p) = -V F (/>f(p) J (pF(pi)S- 1 (PuT)A(p 1 )dT l -W I( p I (p) J ( p I (p 1 )S- 1 (p 1 ;T)A(p 1 )dT l 

- J Y(p,p 1 )£- 1 (p 1 ;T)A(p 1 )dr 1 . (12) 

This equation is conveniently rewritten as 

A(p) = B fX f(p) + B lX i(p) , (13) 



with 



B 



F = -V F J Mp)^~ 1 (p;T)A(p)dT, 

-Wj j ( j> I {p)£-\p;T)A{p)dr, (14) 



and 



Xf{p) = <f>F{p) - J Y(p,p!)£ 1 (p 1 ;T)xF(pi)dr 1 , 
Xi(p) = hip) ~ J Y(p,p 1 )£- 1 (p 1] T)xi(pi)dT l . (15) 
Referring to the relations (11), we observe that 

Xi{Pi) = Xf(pf) = 1 , Xi(Pf) = 0, Xf(pi) = <J>f(pi) = V(pi,Pf)/V(pf,Pf) , 

(16) 

because the block Y is zero when either of its arguments lies on a Fermi surface. By 
this same property, it is permissible, inside the quantity £~ x appearing in the integral 
equations (15), to replace the superfluid quasiparticle energy E(p 1 ) by |£(pi)| and 
the temperature factor tanh [E(pi)/2T] by unity. Because the energy gaps involved 
are generally quite tiny compared to the Fermi energy, these replacements are valid 
to a superb approximation. We are left with the linear integral equations 

Xf(p) = </>f(p) ~ J Y (p, Pl ) 2 ^ 1 xf (pi ) dri , 

Xi{p) = <t>i(p) ~ Jy( P , Pi ) ^L-y XjfyJdT! , (17) 



for the two shape functions needed to construct the gap function A(p) using Eq. (13). 
Since there remains no trace of the temperature T in Eqs. (17), we are free to regard 
the solutions Xiip) an d Xf(p) as T-independent quantities. 
Appealing to the properties (16), Eq. (13) yields 

A F = A(p F )=B F , Aj = A(pj) = B I + B F <f> F {p I ). (18) 

Inserting the decomposition (13) into Eqs. (14), we arrive at a system of two equations 

B F = — VpLppBp — VpLpjBj , 

S 7 = -W I L IF B F - WiL n Bi , (19) 
for determination of the amplitudes Bp and Bj entering Eq. (13), where 

L FF = J <j)F(p)£~ 1 (p;T)xF(p)dT , L n = J (p I (p)£~ 1 (p;T)xi(p)dT , 



jjF = J^FI 



L FT — 



J Mp)£~ 1 (p;T)x F (p)dr = J cj>F{p)£~\p;T)xi{p)dT . (20) 

It is helpful to recast the system (19) in an equivalent form 

[1 + V F L FF (T) - VfMpi)Lfi(T)] A f + V f L F i(T)Aj , 

[WjLFiiT) - (1 + W I L II (T))MPi)] A F + [1 + WjLuiT)} Aj = . (21) 

For a solution to exist, the determinant V(T) of (19) or (21) must equal zero for any 
T. Together with either of the two equations (19) [or either of (21)], the dispersion 
relation T>(T) = forms a closed system that permits one to determine all charac- 
teristics of the superfluid system feeding upon the two Fermi surfaces located at pi 
and pf- 

We begin to explore the implications of the formalism we have developed by 
examining the effect of the additional (inner) Fermi surface on the superfluid transi- 
tion temperature T c . Observe that Eqs. (19) [or (21)] become decoupled if Ljp = 0. 
Assume, as a first case (Case I), that both of the interaction parameters Vp and Wj 
are negative, so that Cooper pairing could exist at both Fermi surfaces when they 
are disconnected. The pairing effect is naturally stronger at the "main" or outer 
surface due to a greater density of states. From the two solutions of the problem as 
stated, we therefore select A(p) = Afxf(p) with A/ = (pF(pi)Ap, implying that 
the individual critical temperatures Tf and T/ satisfy Tf > T/ . It should be noted 
that in spite of this inequality, the magnitude of the ratio Aj/Ap ~ cppiPi) is not 
necessarily less than unity (see below). 

Working in the vicinity of the transition temperature T c , standard arguments 
and manipulations in the spirit of BCS theory reveal the behaviors 

L FF (T - T c ) - N F (0) {(1 + gj F )(L + a x r) - a 2 [D% + gj F Dj] } , 
L n (T -> T c ) -> JV 7 (0) [L + ai r - a 2 Dj] , 

L IF (T -> T c ) -> NtWMpi) [L + «ir - a 2 Dj] , (22) 



in terms of the three dimensionless parameters r = (T c — T)/T c , Dp = Ap/T c , and 
Dj = Aj/T c . In Eqs. (22), Np(0) and iVj(O) are respectively the densities of states at 
the outer and inner Fermi surfaces, gj F = 4> f (pi)Ni(0)/Nf(0) is an effective coupling 
constant, and L = \n(s F /nT c ) + C measures the transition temperature, where e F is 
the free Fermi energy and C is Euler's constant (0.577). Certain constants entering 
the derivation of the limiting behaviors of L//(T), L[p(T), and Lff(T) have no 
material role in our arguments; they merely effect a renormalization of the critical 
temperatures T F and T/ and may be thus be suppressed in forming Eqs. (22). The 
relevant temperature dependences are determined entirely by the ratio a\ja 2 = 
8tt 2 /7C(3). 

After substituting the results (22) into Eqs. (21), we arrive at 
[1 + V F N F (0)(L + a x r - a 2 D 2 F )] D F + VpN^O^piPi)^ + out - a 2 D 2 I )D I = , 

-<t> F {pi)D F + [1 + WjWj(O)] (L + ol\t - a 2 D 2 I )D I = . (23) 

Setting T = T c and evaluating the determinant V(T C ) of this system, we obtain a 
closed formula for the new critical temperature T c in terms of the individual critical 
temperatures T C F and T C J for the uncoupled system: 

(L-l I )(L-l F )-gj F l I L = 0. (24) 

The constants l F and // entering this relation are defined by l F = \n(s F /7rT^) + C = 
-1/V F N F (0) and h = ln(e° /ttT/) + C = -l/WjiVj(0). The inequality T C F > T C J 
clearly implies // > l F . 

The situation at small coupling, gj F « 1, is especially transparent. In this case 
the quantity L (which measures T c ) is not much different from l F (which measures 
T F ), allowing us to replace L by l F in the last term of the determinantal condition 
(24). The solution of Eq. (24) is then 



T Ii + If , 

L+ = ± 



(//-/ 



F 



\'2 



+ g^F^i^F 



(25) 



This result reminds us of the familiar textbook solution of the two-level problem. In 
that problem, the two energy levels repel each other when the off-diagonal interaction 
is switched on. The lower level moves downward and the upper level moves upward. 
In the current problem, the greater logarithm (in this case, L+) increases while the 
smaller logarithm (L_) decreases. In particular, 

L_-l F c- g 2 IF J^-. (26) 
// — ip 

In the case being considered, li and If are both positive, with Zj > lp- We may 
therefore conclude that the emergence of the second Fermi surface increases the 
critical temperature T c relative to T C F . 

The picture changes nontrivially when we consider the more interesting case 
(Case II) in which pairing is absent at the new Fermi surface when the two surfaces 
are disconnected, but is still present at the original surface. When the coupling is 



turned on, a pairing gap is found to exist on the new Fermi surface as well as the old 
one, a feature that is directly evident from either of the equations (21). The result 
(25) continues to apply (again assuming small coupling). However, in contrast to 
Case I, the value of // becomes negative while / f remains positive. Consequently, the 
single acceptable value of L derived from Eq. (25) increases relative to If, implying 
a decrease of T c with respect to Tf . This behavior should not be surprising: the 
value of the pairing gap depends on the shape of the single-particle spectrum, and 
if the spectrum becomes flatter in a region where the interaction is repulsive, there 
must be a suppression of the gap value and a concomitant suppression of T c . We 
must stress that the situation is now quite different from that of perturbation theory, 
where the gap increases independently of the sign of the perturbating interaction. The 
distinctive behavior we have described is indicative of a failure of perturbation theory 
in Case II. We should also point out the close resemblance between the predicted 
behavior and the proximity effect observed in junctions between a superconductor 
and a normal metal: the superconductor tends to induce superconditivity on the 
normal side of the junction, while suppressing its strength on the superconducting 
side. 

Here we have confined our attention to the effect of the bubble rearrangement 
on the superfluid transition temperature T c , considering the possible scenarios for 
pairing when there are two concentric Fermi surfaces. Results are also available [18] 
for the modifications produced in the specific-heat discontinuity at T c and for the 
relation between T c and the energy gap Ap at T = 0. 

4. CONCLUSIONS 

We have studied the rearrangement of single-particle degrees of freedom that 
precedes the onset of pion condensation, and we have found that this rearrangement 
may express itself in the emergence of a bubble in the quasiparticle momentum dis- 
tribution, i.e., the formation of a new Fermi sea with a spherical hole in the middle. 
This is in fact one of the scenarios considered some two decades ago by de Llano, 
Plastino, and their collaborators. The formalism we have developed and the results 
we have obtained can be applied more widely in the theory of strongly correlated 
Fermi systems. In this spirit, it will be of special interest to re-examine the case 
of superfluid 3 He, which offers a realistic example of a Fermi liquid existing near an 
antiferromagnetic phase transition. In our view, rearrangements of the single-particle 
degrees of freedom arising from momentum-dependent effective interactions have a 
generic character. Such rearrangements include not only bubble configurations, but 
also the phenomenon of fermion condensation [8-12]. We call attention especially to 
Ref. [11], where the competition between bubble rearrangement and fermion conden- 
sation is studied under variation of the temperature. 
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